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Abstract 

The electrical resistivity including the Kondo resistivity increase at low tem- 
perature is calculated for thin films of dilute magnetic alloys. Assuming that 
in the non-magnetic host the spin-orbit interaction is strong like in Au and 
Cu, the magnetic impurities have a surface anisotropy calculated in part I. 
That anisotropy hinders the motion of the spin. Including that anisotropy 
the effective electron-impurity coupling is calculated by using the second or- 
der renormalization group equations. The amplitude of the Kondo resistivity 
contribution is reduced as the position of the impurity approaches the surface 
but the increase occurs approximately at the bulk Kondo temperature. Dif- 
ferent proximity effects observed by Giordano are also explained qualitatively 
where the films of magnetic alloys are covered by pure second films with dif- 
ferent mean free path. The theory explains the experimental results in those 
cases where a considerable amount of impurities is at the surface inside the 
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ballistic region. 
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I. INTRODUCTION 



In the previous paper |^ (Part I) following Ref's and we have calculated the 
magnetic anisotropy for a magnetic impurity embodied into a non-magnetic host (e.g. Au, 
Cu) with large spin-orbit interaction for the conduction electrons on the sites of the host 
atoms. The magnetic anisotropy is developed due to the exchange interaction between 
the magnetic impurities and the conduction electrons. As the scattering due to impurities 
involves different angular momentum channels (— / < m < I), therefore, the scattering 
depends on the directions of the conduction electrons before and after the scattering. In 
that way the scattering on the impurity itself depends on which host atoms the electrons are 
scattered by the spin-orbit interaction. Due to that dependence the anisotropy is determined 
by the positions of host atoms around the impurity and it is the larger the stronger the 
asymmetry around the impurity. As that information is carried by the momenta, therefore, 
that is restricted to the atoms in the range of the elastic mean free path l^i- Thus the 
anisotropy can be developed only if the impurity is inside surface area of the thickness of 
the elastic mean free path (ballistic region of the surface). If the surface in that region is 
plane like then the anisotropy energy is (see Eq. (1) of Part I [0) 



where S is the spin operator of the impurity, n is the normal direction of the experienced 
surface element, and Kd is the anisotropy constant depending apart from the oscillatory 
part like l/rf on the distance d, measured from the surface. Of course, if the ballistic region 
contains more sophisticated part of the surface, then the determination of the direction and 
amplitude of the anisotropy is a more complex task. The oscillatory part decays faster than 
1/d approaching the bulk part of the sample (see Eq. (B20b) in Part I). 

In the surface area the anisotropy energy leads to different splitting schema shown in Fig. 
ID, depending whether the spin S is integer or half-integer. In this way the spin very nearby 
the surface is freezing into a singlet or doublet considering the integer and half-integer case. 




(1) 
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respectively. Thus at low enough temperature the spin shows no, or restricted dynamics. It 
is important to point out, that the states m = ±1/2 for S = 5/2 cannot be replaced by a 
doublet of 5 = 1/2. The spins are rather squeezed into planar states as shown in Fig. ^ 

Assuming that Kd ^ T, the integer spin does not contributes to the resistivity in contrary 
with the case of the half-integer spin where the two lowest states contributes to the resistivity. 
The spin S = 1/2 is, however, not affected by anisotropy. Thus different behaviors of 
the electrical resistivity can be expected depending on the value of the spin. Considering 
impurities nearby the surface inside the ballistic region more and more orbitals become 
active as the impurity positions approaching the bulk (see Fig. ^). 

Similar structure appears in the temperature dependence of the resistivity. Cooling down 
the sample, at the beginning almost all the spins are free. At further cooling more and more 
spin states are frozen in thus in case S = 2 aX T < Ti three orbitals and T < T2 a single 
orbital is populated (Ti : T2 = 4 : 1) while in case of S' = 5/2 each state is double degenerate 
(Ti : T2 = 25 : 9). 

Assuming that in the region of the Kondo temperature the occupations of the different 
states are varying by considerable amount depending on the positions of the impurities in 
the surface region, then lowering the temperature less and less impurities can further develop 
the Kondo state, thus less and less impurities can contribute to further increase of the Kondo 
resistivity. As it has already been pointed out the reduction in the contribution to resistivity 
somewhat less pronounced for half-integer spin than for integer spin (see the discussion of 
Fig. The contribution to the Kondo resistivity can be schematically plotted for impurities 
with different distances measured from the surface < d2 < di (see Fig. |). 

The phenomena described above are very similar to the Kondo effect in the presence of 
crystalline field at the impurity site when lowering the temperature different crystal field 
states are frozen out, but those fields are identical for all impurities of the same kind. Such 
calculations have been performed e.g. by Shin-ichi Kashiba et al. 

The reduction in the averaged Kondo resistivity is sensitive on the size of the sample e.g. 
film thickness or diameter of the wire etc. as the ratio of the surface influenced impurities to 
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the total number of the impurities goes to zero as the size is further increased. The role of 
the surface on the impurities can be reduced by depositing another pure film on the surface 
of the sample. The effect can be also influenced by changing the elastic mean free path in 
the samples or in the deposited films. The effects of those kinds will be summarized and 
discussed at the end of the paper (see Sec. |Vli] ) making use the qualitative results obtained 
for the resistivity in Sec. ^ and ^ and all of the references can be found there. 

For the actual calculation of the resistivity first the temperature dependence of the ef- 
fective exchange coupling must be calculated. For that, except the very low temperature 
region, the second order multiplicative renormalization group transformations (two-loops ap- 
proximation) can be used which give smooth behavior at the Kondo temperature in contrast 
to the first order scaling (one-loop approximation) which results in an artificial divergence 
at the Kondo temperature. Those methods are generalized by taking into account the sur- 
face anisotropy terms which occur as a low energy cutoff of the logarithmic integrals in the 
calculation of certain diagrams. The different diagrams depending on their spin labels have 
different infrared, low energy cutoff due to the anisotropy. These calculations are in close 
analogy to those with crystalline splitting. The next step is using these effective couplings 
to calculate the electrical resistivity by solving the Boltzmann equation and an average over 
the impurity positions is taken also. 

Finally, to fit the calculated resistivity at low temperature, an effective surface layer 
thickness A can be introduced, by assuming that inside of that surface region there is no 
Kondo effect, and outside the Kondo anomaly is fully developed. The experimental data 
are compared both with that phenomenological description and the original calculations and 
they are giving equally excellent fit. 

The paper is organized as follows. In Sec. ^ the general scheme of the multiplicative 
renormalization group (MRG) is presented for the Hamiltonian with the anisotropy term. 
The scaling equations are presented in Sec. |T| which are solved in Sec. |rV|. The electrical 
resistivity contribution is calculated in Sec. |V| for the dilute limit and Kondo resistivity in 



thin films is given in Sec. 0. Sec. |VI1| is devoted to the experimental results, theoretical 



interpretation of the results and prediction. A general discussion is contained in Sec. [VIII] . 
The Appendix contains the actual calculation of the diagrams which are needed in Sec. fT\\ . 
Throughout the paper the h = ks = I units were used. 

II. THE HAMILTONIAN AND THE GENERAL FORM OF THE MRG 

TRANSFORMATION 

The Kondo Hamiltonian in the presence of the anisotropy is 

H = Y^Sk O'laO'ka + Ha 

+ E JmM'SmM' {(^\a<^<T(T'0.k'<T')-, (2) 

M,M' 

where a}.^ {dka) creates (annihilates) a conduction electron with momentum /c, spin a and 
energy Ek measured from the Fermi level. The conduction electron band is taken with 
constant energy density po for one spin direction, with a sharp and symmetric bandwidth 
cutoff D. a stands for the Pauli matrices, Jmm'^s are the effective Kondo couplings and Ha 
is given by Eq. For the impurity spin the Abrikosov' pseudofermion representation 
was used 

S = h^MS MM'^M' 1 (3) 

where the projections of the z component of the impurity spin are described by an auxiliary 
fermionic field 6m (M = —S. . . S). Choosing the quantization axis parallel to n, with this 
substitution the Hamiltonian Eq. (^) become 

H = Yjk a-l^aka + + KM)h\jhM 

k,a M 

+ E 'JMM'{b^M^MM'bM'){al^(^o-cT'ak'a'), (4) 
k,k',a-,cr' 
M,M' 

where the chemical potential A oo was introduced to project out the physical pseud- 
ofermion subspace J2b\ibM = 1 and the notation Km = KM'^ was introduced for the MRG 

M 

calculation. 
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The conduction electron and pseudofermion Green's functions are 



LJ — Ek — 



and 



Gmm'[uj) = (6) 

a; — A — Lm(u;) 

where u = u — Km- and Ej\/(ct;) are the self-energies for the conduction electrons and the 
pseudofermions, respectively. They are diagonal in the adequate spin quantum numbers, 
because of that the whole Hamiltonian is symmetric under rotation around the z axis. The 
vertex function is denoted by T k^M,k'a-'M' {^1,^2, ^3,^4) ■ 

The multiplicative renormalization group transformation can be written as ^ 

Gk(7,k'a'{^^, JlVIM'iD), Km{D), D) = Z (^-jj^Gkc7,k'a'i^, JmM'^ ^o) i'^^) 

gMM'{C0,JMM,{D),K^,{D),D) = ZA/(^)6^AfM'(^,4/Af"^°f ^0) (7b) 

r T (n\ K (m n\ - ^'"^M,k'a'M'iuji, J°-^^-^,, Kl^, Dp) 

where Z{^) and Zm{^) are the renormalization factors for the electrons and pseud- 



ofermions, respectively. Introducing x = ln(-^) as scaling parameter, the Callan-Symantzik 



MRG equations are 

^ dG ^ ^ dG ^ dG 
-vG+^+Y. I3mm'^. + E 7M^T^ = (8a) 

•^•^ MM' CI J MM' M CIJ\M 

- m.QM + ^ + E ^MM'Qj—^ + E = (8b) 

AfAf A/A/' M M 

(^ + + 2 w)rA/A/' + + E I^mm'-qJ— + E 7m g^, = 0, (8c) 

MM' MM' M 

where 



d\nZ , , 

r, = — (9a) 

VM = (9b) 

= (9c) 

7M = ^ (9d) 

and for the sake of simplicity the electron and pseudofermion Green's function, and the 
vertex function were denoted by G, Qm, and Tmm', respectively. The initial values for the 
renormalization factors and couplings are Z = Zm = ^, Km = K^m — KM'^, J mm' = <^ for 
each M, M', at D = Do- 
Using the definition of self energies in Eq. @ and (P) the first two equations can be 
rewritten as 

- [uj-Ek- Se)?7 + ^— + 2^ Pmm'^ h 2^ IM^TT- = (10a) 

MM' CI J MM' M OJ\m 
-{U-X- J:m)VM + -TT- + L /^MAf'^TT + L + -T" = (10b) 

which form is more comfortable for calculating the MRG equations. 



III. CONSTRUCTION OF THE MRG EQUATIONS 

To construct the MRG equations the perturbation theory was applied, that is the Hamil- 
tonian was divided into a non-interacting and an interaction part with small parameters 

<^AfA/''s. 

For which the electron self energy contains a closed pseudofermion loop. Eg tends to zero 
as A ^ oo. Thus in the thermodynamical limit for a single impurity from Eq. (|10a| ) 7^ = 
and Z = 1. 
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Turning to the other two equations Eq. (|10b| ), (|8^) they were solved in next to leading 
logarithmic approximation where the MRG equations are 

(i-AR? = ^ + 7l? (11) 

and 

-Q^ r <^aa''^MM'PMM' — U, [LZa) 



/I (2) 1 (2) X J- c- I ^'^^Sv/,<7'M' 
[7:Vm ~r 7:VM')--'MM'(^aa'^MM' -\ 



dx 

MM' f)l - - 



+ H ^i^Af' P)^^'' ^ + (^aa'SlMM'^MM' " (12b) 
MM' "-^MM' 



where rf^j , 7!^'' and I3^mm' proportional to the second, 13'^mm' the third power of 
the Jmm'^s, respectively. The whole next to leading logarithmic /3-function is Pmm' = 
Pmm' ^ Pmm' ■ 

Thus to construct the next to leading logarithmic scaling equations we have to calculate 
the second {S'^o-Mo-'m') third {^'^aMa-'M') order vertex corrections, and the second order 
self energy correction (5E\/) for the impurity spin. 

These corrections were calculated by applying the thermodynamical Green's function 
technique and analytical continuation [0. Assuming scaling for the vertex function Tmm' 
only one energy variable was kept [§], thus ui = u, UJ2 = Km, uj^ = u + Km — Km' and 
UJ4 = Km', where ui, uj2 the incoming, tus, UJ4 the outgoing electron and pseudofermion 
energies, respectively. 

The second and third order vertex diagrams, and the second order correction to the 
self-energy for the impurity spin are shown in Fig. |^ and Fig. |^, respectively. The detailed 
calculation of these diagrams is carried out in the Appendix. Collecting the whole second and 
third order vertex corrections together, they and the self-energy correction for the impurity 
spin were substituted into the Eq's ( pi]) and (|1^). In Eq. (]T2[b) the contributions of the 
third order parquet-type diagrams depicted in Fig. ^ (b) and (d) were canceled out with the 
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?(2) '^^aM.a'M' 



terms J2 P\-Jic„ aj'^".^' as the leading logarithmic scaling equations are equivalent to the 



MM' 



MM' 



summing up of the parquet diagrams. The divergences at finite T were canceled out, too. 
Thus only the diagram Fig. || (c) contributes to Eq. (|l2|b). 

Introducing the dimensionless couplings Jmm' = PqJmm' the next to leading logarithmic 
scaling equations are 

riM = ^^^^^ = q^(S, M)jM,Af-ljM-l,M0M,M-l(-D) 

ax 

+ p\S, M)jm,m+iJm+i,m^m,m+i (D) + ^M^jli^M (13) 

7m = ^^^^^ = {Km-1 - KM)q^{S,M)jM,M-ljM-l,M'^M,M-l{D) 

+ {Km+1 - Km)p^{S, M)jM,M+dM+l,M^M,M+liD) (1^) 

Po(3m,M+1 = '^^^^^^ = -{MjMM - (M + l)jM+l,A/+l)jM,M+l(l + 0M,M+l(-D)) 
+ q^{S, M)jM,M-ljM-lMjM,M+lQM,M-li^) (^©Af.M+l - 

,M+ljM+l,Af+2jM+2,A//+10A/+l,Af+2(-D) (^M,M+i 

- ^^(5*, M)j^ j^,^^;^jM+l,A/0A/,Af+l(-D) — M'^ilj^j^jM^M+l 

- (M + l)^jM+l,A/+liM,M+l + 2M(M + l)jM,AfjM,Af+ljM+l,Af+l (15) 
Po(3m,M-1 = = -((M - l)jM-l,A//-l - MjM,M)jM,M-li^ + 0Af,A/-l(^)) 

ax 

+ p^{S, M)jM,M+dM+l,MjM,M-l^M,M+liD) (^0Af,Af-l - 

+ q^{S,M - l)jM,A/-ljA/-l,A//-2jA/-2,A/-10Af-l,A/-2(^)(^0A/,A/-l(^) " ^) 

- ^^(5", M)jlj j^^j_jM-l,M^M,M-liD) - 

- (M - 1)^jX^_i_^_iJm,A/-1 + 2M(M - l)jM,A^jM,Af-ljA/-l,A/-l (16) 

and for M 7^ 

^ _ djM,M 

PoPm,m — — J — 
ax 

q^{S, M)QM,M-liD)jM,M-dM-l,M - P^{S, M)&M,M+liD)jM,M+dM+l,M 



M 
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+ (fiS, M) — — — jM,M-ljM-l,A/-ljAf-l,M0M,M-l(-D) 



+ p {S,M) — — — jM,M+ljM+l,Af+ljA/+l,M0Af,M+l 

- q^iS, M)iM,M-lhl'l,Mhl,M^M,M-l{D) 

- p^iS, M)jM,M+ljM+l,MjM,M^M,M+l{.D) (17) 



where Qmm'{,D) ensures that D > J {Km — Km>Y + T"^ with the definition 



_ 1 ilD> J{Km-Km'Y + T'^ 

Qmm'{D) = { \ (18) 

iiD<J{KM-KM'Y + T\ 



The definition of ^(5*, M) and q{S, M) are given in Eq. ([A7| ). It must be stressed that these 
scaling equations are vahd for T < D. 

IV. SOLUTION OF THE SCALING EQUATIONS 

It can be seen from Eq. ([T^), that the corrections to the bare K^^ are proportional to 
the second power of Jmm' 's in the leading order. In the other equations Km appears only 
in the arguments of the G functions which are multiplied by the second or third power of 
Jmm''s. Thus to keep in the approximation consistent Km = Km = KM"^ can be taken 
in the arguments of the G functions. After this consideration the solution of the equations 
becomes simpler, because the equations for [3 mm' are not coupled. Exploiting the symmetries 
of the scaling equations 

iM,M' — jM',M (19a) 
jM,M' = j-M-M' (19b) 

must hold. 

Thus the equations which have to be solved are 

djlMM+l ^ -(MjM,M - (M + l)jM+l,A'/+l)jM,A^+l(l + Gm,A'/+i(-D)) 

ax 

11 



+ q (5',M)j^^jv/-liM,M+10M,M-l(^)(^0M,M+l(^) - 2 
+ p'^iS, M + l)jAf,M+ljM+i,M+20A/+l,A/+2(£') (oMM+li^) " ^) 

- (M + l)^jM+i,M+liM,M+l + 2M(M + l)jM,MjM,M+ljM+l,M+l 



(20) 



and for M 7^ 



g2(5,M)eM,M-i(/^)jl,M-i 



djM,M _ 1 

-/(^,M)eM,Af+l(/^)jM,M+l 

M- 1 



Jm, Af - iJM - 1 , Af - 1 Af , A/ - 1 



M 

^ + 1-2 ■ /n^ 

— j^Jm,m+iJm+i,m+i'^m,m+i{^) 



- q^iS, M)jlj j^_jM,M^M,M-liD) 

- p^{S, M)jlj j,j_^^jMM'^M,M+l{D) 



(21) 



These equations were solved numerically for different initial couplings jo = jMM'{D = 
Dq). The results for jo = 0.1 and jo = 0.0435 at K/T = 10 are shown in Fig. |^ and Fig. || for 
S' = 2 and S = 5/2, respectively. The initial bandwidth cutoff was chosen as Dq = 10^ K. 



V. RESISTIVITY 

The Kondo resistivity was calculated by solving the Boltzmann equation in the presence 
of the spin-orbit induced anisotropy, using the value of running couplings {Jmm') calculated 
in the preceding section, at D = T. 

Taking the usual form for the electron distribution function /(e^) in the presence of 
the electric field E as 

f{e,) = foisk) - (kE)$(£,)^:M£^, (22) 

OEk 

the Kondo contribution to the resistivity is 
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' ^ (2m)3/^/&..f$(..)(-^) (23) 



where the function $ is determined by the Boltzmann equation 



£^(UE) + (I) 

m OEk \ot J 



coll. 



0. (24) 



The colhsion term (§f )cou can be expressed in terms of transition probabihties as 

-W'(k,ff^k',ff')/(k)(l-/(k'))} (25) 

where e.g. (k', a' k, a) represents the transition probabihty from a k', a' state to a 
k, (7, c is the impurity concentration, and V is the volume. 
Turning to our case, these probabihties can be calculated as 

Vr(k, a k', a') = J2 PMw{k, a,M ^ k', a', M') (26) 

MM' 

where 

w;(k, a, M k', a', M') = 27r|Tk,.,M^k',a',Aff5(£fc - + KM^ - KM'^), (27) 
and pm = e-^^*^7z, /5 = l/T, Z = E e"^^^'. 

A/ 

The scattering amplitude in Eq. ( [27| ) is expressed in terms of the renormalized couplings 

Jmm'{x) as 

T]i^cT,M^'k',cT',M' ~ JmM'{x = \'il-^)(Tcrcr' S MM' (28) 

where the dependence on the direction of the momenta k and k' is ignored and that makes 
the Boltzmann equation solvable in a simple form. The k-dependence may result as some 
numerical factors in the final expression, but in the main features of the temperature de- 
pendence those are not playing an important role. 

Substituting these assumptions into Eq. (p3), changing the sum to 

k' 

J p{ek')dek' J ~ Po / de J using the properties of the spin algebra for cr and S 

-Do 

and the "detailed balance" principle (( ^'^'gf^'°^ )coii ~ obtain after linearization in E 
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where 



Fi^k) = Y.PM\M'^jljj^{x = In 

M *■ 



Do. 
T ' 



+ (SiS + 1) - M(M + 1))j1,m+i(^ = In ^) 

■[1 " (1 - e-(^*^+^)^/^)/o(.. - (2M + (30) 



where we introduced the dimensionless couphng constants Jmm' = PqJmm'- Inserting Eq. 
(p9|) into the Boltzmann Eq. ( p^ we obtain for $ 

-2ttc 



Hsk) = -- —Fisk) 
m I po 

and for the Kondo resistivity 



(31) 



' ' '■*f-^)^'-(-) (32) 



PKondo P^°^ J \ de 
where the usual assumption 

dfo\ X _ 3/2 f ,J df( 



/^-'1-f)^-(^)-^'7'^(-f)^-'(^) (33) 



was taken into account and the constant = It^^^ was introduced. 



In the case of = 0, Eq. ( p2D reproduces the bulk Kondo resistivity. 

The resistivity was calculated by evaluating the occuring integral in Eq. (^) numerically 
for different K values which are in the regime discussed in Part I The resistivity as the 
function of the temperature is shown in Fig. ^ (a) and (b) for S = 2 and S = 5/2, respectively. 
The plots are similar to the experimental ones (see Sec. |Vli| ). 

The effect of the anisotropy on the Kondo temperature defined by the largest slope in 
the resistivity can be examined by looking at the derivative of the calculated resistivity vs. 
temperature. We can see from Fig. ^ (a) and (b) that the Kondo temperature defined in 
that way is only slightly affected by the anisotropy in those cases where the Kondo effect is 
pronounced (e.g. K < 0.5 K in Fig. The effect of the anisotropy becomes dominant for 
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larger strength of K. The temperature dependence of the resistivity has a maxima depending 
on the strength of K and the spin-flip contribution is freezing out gradually. That behavior 
is very different for integer and half-integer spins for large anisotropy. For integer spins 
the impurity contribution tends to zero in a way which is very sensitive on the strength 
of the anisotropy. In the case of half-integer spin the impurity resistivity approaching, 
however, a finite value at zero temperature which is independent of K. The resistivity there 
is determined by the dynamics of the two lowest energy levels shown in Fig. ^. In those cases 
the Kondo effect is also essentially reduced due to the smallness of the spin-flip amplitudes, 
but still presents. (The Kondo contribution in the lowest order is proportional to Scr which 
gives a small amplitude for M = ±1/2.) 

It is important to emphasize, that in case S = 1/2 the anisotropy is loosing its meaning. 

In a real system the anisotropy strength K has a distribution thus the formation of the 
resistivity maximum at finite value of temperature cannot be expected at least above the 
Kondo temperature. The calculation is anyhow not reliable for T < T^- 



VI. KONDO RESISTIVITY IN THIN FILMS 

To get some information about the case of thin films a simple assumption is made that 
the two surfaces contribute to the anisotropy constant K in an additive way. The anisotropy 
factor for a sample with thickness t and in a distance d measured from one of the surface is 

rv rv 

K{d, t) = Kd + Kt-d = ^ + j—^ (34) 

where the coefficient a is estimated in Part I [|l| (see Eq. (32)). The appropriate calculation 
of the resistivity including the elastic impurity scattering with mean free path /ei is a very 
difficult task for a film for an arbitrary ratio of lc\/t and value of K. In order to avoid those 
difficulties we are making use of the fact, that the magnetic exchange (Kondo) contribution to 
the resistivity pKondo is smaller by a factor 10~^ than the residual normal impurity resistivity 
Pnor (p = Pnor + PKondo), thus an cxpausiou in the Kondo contribution is appropriate. The 
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calculation can be carried out in two limits (i) t < Id and (ii) t ^ Id- It will be shown 
that the final expression does not depend on which limit is considered. In the case (i) the 
electrical resistivity contains the average value of the inverse electron life time. Denoting the 
resistivity at temperature T for a given value of K by p{K,T), the average over the value 
of K{d, t) is 

1 } 

p{t,T) = -J p{K{x,t),T)dx. (35) 



On the other hand, in case (ii) the sample can be considered as a set of parallel resistors of 
equal size, where each resistor represents a stripe in the sample with a constant K. In that 
case the conductances are additive, thus 

Pit,T) = , 1 „ ^,r^. — TT^ (36) 

i 

where N is the number of the resistors (stripes) labeled by i and 0", represents the Kondo 
conductivity of stripe i placed in distance Xj. In the actual case only the first stripes depend 
on the surface anisotropy. The Kondo conductivity is defined by the Kondo resistivity given 
by Eq. (H) as 



— 0"nor + CTKondo — (Pnor + PKondo) ^ (37) 

where cxKondo ~ — fisgsdo ^ q_ ^jj-^g expansion gives the final expression 

Pnor 

1 } 

p{t, T) = Pnor + J J PKondo{K{x, t),T)dx. (38) 



That expression valid in the limit pKondo ^ Pnor gives back exactly the expression in Eq. 



In the numerical calculation the integral in Eq. (RH) or (p3) is replaced by a weighted 



sum with appropriate intervals. Introducing the new integration variable d/a, the calculated 



Kondo resistivity depends only ont/a which is shown in Fig. |10| (a) and (b) for 5 = 2 and 
S = 5/2, respectively. Fitting the calculated Kondo resistivity for temperatures T 3> 
Tk (T = 2 — 4 K) by the function pKondo/p^^^ = — -BcaiclnT as it has been done in the 
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experimental works (see Sec. [VI I|) , the behavior of i?caic was examined as a function of t/a 
which can be seen in Fig. |Tl|. To compare this calculated dependence of the coefficient B on 
the thickness to the experimental data, they were fitted by the function B{t) = p^^^Bca\cit/a) 
as it is shown in Fig. ^ The fitted value of a is a = 247. 7A K which is in agreement with 
the prediction given in Part I by Eq. (32) (see Sec. |Vli| ). The fit is not too sensitive to 
small changes (< 5%) in a. 

If the sample is not thin, then the above results can be phenomenologically described in 
the framework of a simple model where the impurities in the region of the surface do not 
contribute to the Kondo resistivity and outside that region they are not affected. In this 
way the effective suppression length A can be introduced and then the average resistivity at 
low temperature pt e.g. for a thickness t 

t-2X 

Pt = Pt=oo — - — • (39) 

According to this semi-phenomenological formula B{t) = Boo{l — 2X/t) which was fitted 
to the experimental data. This can be seen also in Fig. where the fitted value of the 
effective suppression layer parameter is A = 207. 5 A. 

The effect of the mean free path in the ballistic region can be demonstrated directly 
by taking into account the effect of the mean free path in the anisotropy constant. We 
calculated the change of the electrical resistivity for a thin film with thickness L = 600A 
with anisotropy arising only at one of the surfaces in the forms 

K = A-^ (40a) 
a 

K = A- t-'^l^ (40b) 
a 

where ^ is the elastic mean free path (e.g. ^ = lOOA), d is the distance measured from the 
surface with anisotropy of strength A, and the exponential decay is due to the mean free 
path. The electrical resistivity is calculated for S* = 2 at T = 2Tk just above the Kondo 
temperature T^, as a function of the strength A of the anisotropy for two cases without and 
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with exponential factor (see Eq. (0)). Increasing the anisotropy strength A the spins are 
completely frozen in nearby the surface, but that region is limited by the finite mean free 
path. Fig. |T3] clearly demonstrates that the strength of the anisotropy and the size of the 
suppression layer are reduced due to the finite mean free path as it is calculated by taking 
into account the anisotropy only for one of the surfaces. 



VII. COMPARISON WITH EXPERIMENTS 

In the last couple of years a very extensive study of the Kondo effect in thin films and 
wires have been performed. The experimental works were concentrating on determination 
of the effect of reduced dimensions on the Kondo temperature Tk and the amplitude of the 
resistivity anomaly. A detailed critical discussion of the earlier works are given in [|10|. The 



early studies have been performed by Giordano and his collaborators [^,|T^, and by DiTusa, 
Lin, Park, Isaacson and Parpia [|1^]. In order to discuss the effect of uncoupled magnetic 
impurities, first only those experiments are listed which are performed in dilute limit, thus 
e.g. for Au(Fe) alloys the Fe concentration is 30 ppm. These experiments belong to two 
groups depending whether size effect was observed or not. 

Concerning the theory two regions must be distinguished. When the size of the sample 
(e.g. the thickness of the sample) is inside the ballistic region, then obviously the present 
theory must be applied. In the case of thicker samples more care must be paid. There is 
another theory by Martin, Wan and Phillips which is applicable in the opposite limit 
of weak localization, where the disorder- induced depression or enhancement of the Kondo 
effect is predicted depending on the value of the spin flip scattering rate (depression 
is the case where Tk, hT]7^ ^ T). The competition between these theories needs further 
studies. 

In the following the discussion is organized according to different effects. First we discuss 
how the change in the density of states at the surface can influence the Kondo effect, but it 
is ruled out as an explanation of the size effects to be discussed, because it is applicable only 
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on much smaller scale (Sec. [VII A]) . Then the experiments with considerable dependence on 



the size of the samples are discussed and compared with the present theory (Sec. |VllIj| ). 



Finally those experiments are listed where no size effect was observed (Sec. |V11 C[ ) or the 
concentrations of the impurities are in the spin glass region (Sec. |V11 D| ). 

A. Density of states effects 

As it has been discussed in Sec. | of Part I [|l|, the size dependence cannot be expected 
just because the Kondo cloud cannot fully develop in all directions by reducing the size of 
the sample. The only possibility which has been discussed by Zarand is, that nearby 
the surface there is a change in the density of states of conduction electrons by formation of 
a Friedel type oscillation due to the surface. That explanation was ruled out, because those 
changes in the density of states are very much localized in a few atomic distances measured 
from the surface and the smallest sizes in the experiments to be discussed are about 300A. 
That effect may, however, show up in point contact experiments where the contact size is 
smaller by even more than one order of magnitude. Such experiments were performed by 
Yanson and his collaborators [p!^ -[T8| with Mn and Fe impurities in Cu contacts. Zarand and 
Udvardi [[T9ip0| showed that depending on the actual position of the impurity the density of 



states for an essential energy range around the Fermi surface can be enhanced or depressed 
by even 20%, thus p = po + ^P, where |5p/po| < 0.2. Just in order to demonstrate the 
effect an energy independent 5p is assumed and for that case in the expression of the Kondo 
temperature exp[— l/2J(po + 5p)] = exp[— l/2Jpo] exp[(2 Jpo)^"*^^] there is an enhancement 
due to the second factor. Depending on the value of (Jpo)~^ that enhancement can be over 
a factor of 100 for Mn and about 2 — 3 for Fe impurities. The enhancement is the larger 



the smaller the Kondo temperature Tk [0-|T8|,|20| . In the experiments the enhancement 



is the larger the smaller the contact size, thus to have large enhancement most of the 
impurities must be nearby the surface. Similar effect was also seen E^ in point contacts with 



presumable tunneling two-level systems (TLS) where an atom jumps between two positions 
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and the orbital Kondo effect is developed p2| , p3| by coupling the conduction electrons with 
different angular momenta to the TLS. As the typical size of the studied films and wires are 
much larger and such a dominating enhancement of the Kondo temperature has never been 
observed, therefore that explanation can be ruled out. 



B. Experiments with observed size effect 

Giordano and his collaborators (see for a review [|T^]) have performed a series of different 
experiments under different conditions where the size effect was observed but the changes in 
the Kondo temperature were almost negligible. The experiments of different type are listed 
below. 

(i) Dependence on film thickness. 

The film experiments with thickness 265A — I8OOA were performed e.g. with 30 ppm 
Fe in Au, but similar results are obtained also for 100 ppm 0,|l^ . The resistivity was 
fitted by the formula 

p{T)=p-B\n{T) (41) 

where B is an adjustable parameter. It is well known for the Kondo systems that B 
is just not the result of the first non-vanishing third order perturbational result where 
B would be -B ~ J^, but it is the actual slope nearby or somewhat above the Kondo 
temperature (see for example 1|10[)- In the actual experiments the temperature range 
1.8 — 4 K was studied while Tk = 0.3 K. The dependence of that coefficient B on 
thickness was plotted as shown in Fig. |1^. The experimental results are fitted by the 
calculated dependence of B on the thickness with parameter a = 24:7. iK K and by the 
semi-phenomenological formula given by Eq. (^) with the effective suppression layer 
parameter value A = 207.5A in Fig. |T2| 12|J2§. That value of a is in agreement with 
the estimate given in Part I by Eq. (32). There was not any signal for essential 



change in the Kondo temperature |jTO| in agreement with our theoretical result. It is 
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interesting to note that the estimated Kondo coherence length was about 3-lO^A, much 
larger than the thickness of the sample. Similar experiments were performed with wires 
where more geometrical effects are expected, but the results are qualitatively similar 
but not identical. The simple semi-phenomenological formula given by Eq. ( ^Of ) is not 
appropriate in those cases. Qualitatively very similar results are reported in but 
there are quantitative differences very likely due to the sample preparation. 

(ii) Kondo proximity effect. 

A set of experiments p7|j28[] were performed where the film of dilute alloys is covered 



by a second layer of pure metal. The observation was, that in the case of a thin layer 
of dilute alloys with a significant suppression of the Kondo effect the covering by a 



second pure film results in partial recovery of that suppression. In Fig. |1J with the 
suppression layers indicated it is shown that the bilayer structure has a suppression 
layer only on one side of the film of dilute alloys, thus only one half of the suppression 
is expected. In order to verify the importance of the role of the spin-orbit interaction 
in the superimposed layer to complete the neighborhood of the impurity with a uni- 
form spin-orbit coupling, we suggest experiments where the superimposed layer has 
negligible spin-orbit interaction (e.g. Al or Mg). In that case the boundary is changed, 
but the anisotropy should remain. 

(iii) Kondo proximity effect with overlayers with different disorder. 



It has been shown experimentally that the Kondo resistivity suppression in a film of 
dilute alloys covered by a pure film but with different disorder depends on the disorder 
in the overlayer. It was found that the larger the disorder the smaller the recovery 
is. As it is discussed above if the thickness of the overlayer and the mean free path 
lei in it are larger than the thickness of the suppression layer A, then the depression 
takes place only on one side of the film of dilute alloys. On the other hand, if the pure 
overlayers contains disorder, then the electron entering that overlayer cannot bring 
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back information to the magnetic impurity by their momenta, as their momenta is 
changed in the overlayer (the overlayer is not in the balhstic regime). In these cases 
the reduction in the anisotropy is only partially developed as the surroundings of the 
impurity is not perfectly spherical in contrast to the case of overlayer with long mean 
free path. 



C. Experiments without size effect 



In contrast to the measurements discussed in VII B there is a series of experiments by 



Chandrasekhar et al. where the size dependence was not found. The geometries of these 
experiments were different, the thickness of the sample t was kept the same {t = 380A), 
but the width of the stripes w was changed between 380 — lO^A (see Fig. |1^). After the 
correction due to the weak localization effects and due to electron-electron interaction no 
size dependence was claimed. On the basis of the present theory, for samples t <ti w no 
size dependence is expected as the ratio of the volume of the suppression layers to the total 
volume is not changing. Where t ~ w the anisotropy due to the geometry becomes more 
complicated, thus it is hard to make comparison with the present theory. On the other 
hand for t w the experimental points are somewhat falling off from the main averaged 
line, that, of course, may be due to experimental errors. According to the present theory 
the averaged Kondo resistivity for w ^ t had to be smaller than the bulk resistivity, but it 



seems to be not the case 31 



Finally it should be mentioned that no size effect was observed studying Lai_xCex films 
where Ce has S* = 1 /2 for which no surface anisotropy is expected l3^ . 



D. Higher concentration 



There are several experiments with higher impurity concentration. In these cases 

the impurity- impurity interaction mediated by the RKKY interaction competes with the 
Kondo effect. In another set of experiments |0 the thickness of the film were changed 
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in samples made of Cu with 1000 ppm Cr and it was found very similar depression of the 
Kondo effect described above in Sec. |V11B|. The wires with geometries similar to those 



discussed in Sec. |V11 C] but with 2800 ppm impurities do not show dependence on the width 
d IQ, but the overall amplitude is substantially SU-pressed compared to the bulk, which 
was attributed to spin-glass effects. 

VIII. CONCLUSION 

In the present paper the influence of the spin-orbit induced surface anisotropy is studied 
on the Kondo effect in dilute magnetic alloys samples of finite size at least in one dimension. 
That anisotropy splits the energy levels for impurity spin S > 1/2. That anisotropy reduces 
as the bulk part of the sample is approached relatively slowly as 1/d where d is the distance 
of the impurity measured from the surface. That anisotropy occurs for samples of any shape, 
but for those cases further theories should be developed. The range where the anisotropy 
is relevant can be characterized by the suppression length A introduced in Sec. |VI|, which 
is proportional to the strength of the anisotropy but limited by the mean free path of the 
electron as the anisotropy reflects the presence of the surface in the ballistic region nearby 
the impurity. Thus that suppression length cannot exceed a few hundred A in accordance 



with the experiments discussed in Sec. VII B 



That anisotropy hinders the motion of the impurity spin 5 if 5 > 1/2 and the Kondo 
effect is affected in those regions of the samples where the anisotropy is not negligible relative 
to the Kondo temperature Tk- In order to calculate the Kondo resistivity the renormalized 
exchange coupling constants are calculated in Sec. |II| and |^ by using the multiplicative 
renormalization group technique which is applicable only for temperatures T larger than the 
Kondo temperature Tk, thus no detailed prediction can be made outside that region. It can 
be accepted, however, that if the Kondo effect is already reduced in region T > Tk similar 
effect is expected also for T < Tk- The resistivity is calculated by solving the Boltzmann 
equation in Sec. [V| for integer and half-integer spins with different anisotropy strength. Even 
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if the calculated resistivity curves in Fig. ^ (a) and (b) show different characteristic features 
by developing a resistivity maxima at different temperatures and of different amplitudes, 
those features are almost loosen as an average over the strength of the anisotropy is taken 
for T > Tk- The curves calculated for thin films (see Fig. |T0|) show smooth increase of the 
resistivity. More structures could be expected only in those experiments where the impurities 
are in certain distance measured from the surface. If the anisotropy does not dominate the 
complete sample then as the result of the average taken the largest resistivity slope as a 
function of temperature is in the region of the Kondo temperature Tk, and its position 
cannot be shifted too much on the scale of the Kondo temperature Tk- That theoretical 
result is in accordance with the experimental findings (see Sec. [VII B]) . 



The relatively weak sensitivity of the observed region of the largest resistivity slope on 
the size of the samples rules out the density states effects nearby the surface in contrast to 
the point contact experiments (see Sec. [VII A| ). The size dependence associated with the 
large Kondo compensation cloud is not observed in agreement with the Kondo theory where 
such a simple connection is ruled out. 

The calculated Kondo resistivity for thin films was fitted for temperatures T ^ Tk [T = 
2 — 4 K) by the function pKondo/p^^^ = — -BcaiclnT which is compared to the experimental 
data in Fig. |12| and gives excellent agreement. 

The phenomenological theory using the effective suppression length A (see Sec. [V^) works 
remarkable well to interpret qualitatively the experimental data quoted in Sec. p/II| . The fit 
of the experimental data is shown in Fig. |12l 



The different proximity effects described in Sec. [VII B[ can be also well explained by the 
present theory. 

It is important that the role of mean free path (Sec. |VJ see Fig. |1^) reduces the effect of 
the large anisotropy constant, thus for a large range of strong anisotropy the size dependence 
remains in a limited range as far as the elastic mean free paths are in the same order of 
magnitude. In this way the size effect can be comparable for different host materials with 
different large spin-orbit interactions but with comparable elastic mean free path. 
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We have to emphasize, however, that our calculation does not consider the localization 
effects which are present in samples of larger sizes. Such effect have recently been predicted 
by Martin, Wan and Phillips and deserves further detailed studies. 

In addition to those localization effects the theoretical studies must be extended to the 
microscopic calculations of the anisotropy constant. 

Considering further experiments the mean free path effects should be studied. The most 
relevant experiment to directly verify the role of the spin-orbit interaction could be the prox- 
imity experiments where the superimposed layer would be made of another metal without 
spin-orbit interaction as it is discussed in Sec. |VII B| . In those cases the uniform surrounding 
of the impurity would not be developed, thus the anisotropy remains. Furthermore, the 
experiments with impurities in a certain distance measured from the surface would be also 
very instructive. 

Summarizing, the presented theory is able to provide a coherent description of the size 
effects of the Kondo resistivity in thin films, which is not related to the size of the Kondo 
compensation cloud in any sense. 
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APPENDIX: 

Here we calculate the second and third order vertex corrections, and the second order 
self-energy correction for the impurity spin shown in Fig. ^ and Fig. ^, respectively. 
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Carrying out the Matsubara's summation, analytical continuation, changing the integrals 
/ to / p{e)de f using the assumption for p{e) in Section |, the contribution of 

the second order diagrams are 

^ 1 ( ) 



e -uj + Kj^j" - Km 

iF — ' 

fc — Lu -t J\ M' — -IV ».f" 



D 



-D 



e-uj + Km' - Kj^j 

for diagrams corresponding to Fig. |^ (a). The third order diagrams' contributions are 

^ {l—nF{e)de) ^ {1 — nF{e'))de' 



uj - e + Km - Kn J uj - e' + Km - K^' 

— {(y''cr-'a^)fj(j'JMNJNN'JN'M'SM^Sj^^,S'^/M'Pl 

nAe)de ? nAe')de' 



^ e - UJ + Km' - Kn e' - uj + Km' - Kn' 



(A2) 



for diagrams corresponding Fig. ^ (b), 

— Tr{a^a^)a'^^, JmnJnn'Jn'ai'Sm^SIj,m'S'^n/pI 

^ ^ nF{e){l — np{e'))dede' 



-D -D 



{e-e' + Km - Kn) {e - e' + Km' - K^' 



(A3) 



for diagram Fig. |^ (c), and 



{(J^Cr-'cr'')a(7'JMNJNN'JN'M'SMNSNN'^N'M'Po 

np{e){l — np{e'))dede' 



_^ _^ iio-e' + Km - KN')i£ - e' + Km - K^) 

+ {<J^O'^ 0''^)(t(7'JmnJnN'Jn'M'Smn^NN'^%'M'Po 

^ ^ {1 — nF{e))np{e')d6de' 



_^ _^ {e'-uj + Km' - KN')[e' - e + Km - K^) 

{(T^(y^(y^)ufjiJMN'JN'NJNM'S\iNSNN'SN'M'P'o 

^ ^ np{e){l — nF{e'))dede' 



_^ _^ iuj-e' + Km - KN)ie - e' + Km' - Kn' 

+ {cr''a^a'^)cr(T'JMN'JN'NJNM'SMN^NN'^N'M'Po 
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/ / {1 — nF{e))nF{e')d6de' , , 

'_i _y {e'-uj + Km' - KN){e' -e + Km' - K^') ^ ^ 

for diagrams corresponding Fig. ^ (d). 

The second order correction to the self-energy for the impurity spin according to Fig. ^ 

is 

D D 

dr 

e - e' + Km' - Km 



rr I i j\j T Qi Qi -2 f ^ f ^ ' - nF{e))nF{e') 

-Tr{a a^)JMM'JM'MSMM'SM'MPo / de / de — — — . (A5) 

J J e — e~ 00 + K m' — K m 



For the same indices summation must be carried out. 

The spin factors in Eq. (|A1|) , (|A2|) , ([A^), (^) and (|A5|) were calculated by using the 
identities 

{a'a^)aa' = Sij6aa' + isijka^^,, (A6a) 

(dW'')^^/ = Sija'^^, + 5jk(y'aa' - ^ik(^ia' + i^ijkSaa', (A6b) 
introducing the operators in a usual way, and exploiting that their matrix elements are 



StiM' = Pis, M')8mm'+i = ^/S{S + 1)-M'{M' + 1)6m,m'+i (A7a) 
Smm' = QiS, M')6m,m'-i = ^S{S + 1)-M'{M'-1)6m,m'-i- (A7b) 



Turning to the integrals in Eq. (|Al|), (|A2|) , (|A^), (|A4|) and (|A5|) , after changing the 



integration variable in integrals containing (1 — ^^^(e)) or (1 — nF{e')), from e {e') to —e 
(—£:') they were evaluated in logarithmic approximation. The integrals in Eq. ( [Al| ) give 
logarithmic contribution 



/£,„(fl)=/<i. "-^'1 



-D 

In 



e + uj + Km — Km" 
u + Km - Km" 



D 



+ I(uj + Km-Km"), (A8) 



foT D > \uj + Km — Km"\ and 
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D 



-D 



In 



e — u + Km' — Km" 



uj - Km' + K 



D 



+ liuj - Km' + Km" 



ioi D > \ — UJ + Km' — Km" \ ■ 

The integral in Eq. (^^) gives logarithmic contribution 



D D 



^iiM'i^) 



de j de' 



e + e' + lj-X + Km- Km' 

u) — A + Km — Km' 



D 



-D -D 

^ (cj - A + Km - Km') In 

for D > |a} — A + Km — Km' \ ■ 

The integrals in Eq. (|A2|) give logarithmic contribution 



+ liu + Km - K 



UJ + Km - Kn 



In 
•Mn 



D 

u + Km-K^ 



N, 



N' 



D 



+ I(oj + Km - Kn' 



for D > \oj + Km — Kn\ and D > \uj + Km — K^'l and 



^iJ'NN'i-^) — ^AI'Ni-^) -^AI'N'i-^) 

UJ - Km' - Kn 



In 
■Mn 



D 

UJ — Km' — Km' 



D 



+ + Km -Kn 
+ li-uj + Km' - Kn' 



(A9) 



(AlO) 



(All) 



(A12) 



for D > \ — UJ + Km' — Kn \ and D > \ — uj + Km — Kn' \ ■ 

The integral in Eq. (^^) gives logarithmic contribution for M ^ N , M' ^ N' and 
Km' - Kn' -Km + Kn^O 

n 

np{E)np{E') 



D D 



^ilNM'N'i^) — J '^^ J '^^ ] 



-D -D 

Km — Kn 



e + e' + Km- Kn){e + e' + Km' - Kn' 
Km — Kn 



Km' - Kn' - Km + Kn 

Km' — Kn' 

Km' — Kn' — Km + Kn 



In 
In 



D 

Km' — Kn' 
D 



(A13) 



28 



for D > \Km -Kn\,D> \Km' - Kn'\. For Km' - Kn' = Km - Kn ^ 



Km - K^ 



D 



for D > \Km - Km\. And 



Tie) 

'■MMM 



'M'iD) 



In 



(AM) 



(A15) 



The integrals in Eq. ( |A^ ) give logarithmic contribution 
Imn'mn{D) = I de I de'- 

-D -D 



{e' + u + Km- KN'){e + e' + Km - Kn) 



u + Km — Km/ 



D 



+ I(uj + Km - KN')\n 



(Km - Km? + 



D 



(A16) 



for D > \uj + Km — Kn'\ and D > ^ {Km — K^Y + T^, where T is the temperature, 

D D , . , 



^M'N'MNi^) — J de J de 



nF[e)np[e 



-D -D 

-CO + Km' — Kn- 



to + Km' - KN'){e + e' + Km - Kn] 
+ I{-uj + Km' -KN')\n 



{Km - KnY + T2 



for D > 



D 

uo + Km' - Kn'\ and D > ^{Km - KnY + T\ 

D D 



D 



(A17) 



^\Jnm'N'{^) — j de j de ^ 



nFiejnpie 



-D -D 



lo + Km- KN){e + e' + Km' - Kn' 



'■MNM- 



'N'iD), 



(A18) 



and 



D D 



^ii'NM'N'i-^) — J de J de ^ 



-D -D 



e' -UJ + Km' - KN){e + e' + Km' - Kn') 



^WNM'N'i^)- 



(A19) 



In the estimations above, the function I {a) = J denF{s)^rr^ was introduced which is 

^ ° 

related to the finite T divergences. From the scaling equations the /(a) function is canceled 
out. 

For the sake of handling these contributions more comfortable uo^uo ^ T were set in the 



arguments of logarithms in a way that substituting uo + a {uo + a) with ya? + T^. 
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FIG. 1. The splitting schema due to the anisotropy for (a) S = 2 and (b) S = 5/2. 
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a) b) 
FIG. 2. The squizing of the spin into planar states due to the anisotropy is illustrated for (a) 

S = 2 and (b) 5 = 5/2. The lowest energy states are shown by heavy lines. 
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FIG. 3. The different layers of the impurity positions where more orbitals become active as the 

impurity positions approach the bulk assuming that di,d2 < hi- (a) S = 2, (b) 5 = 5/2. 
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Tk In T 

FIG. 4. The schematic plots of the contribution to the Kondo resistivity for impurities with 
different distances ds < d2 < di measured from the surface. 





a) 



b) 



c) d) 

FIG. 5. The second and third order vertex corrections. The double line represents the impurity 
spin with the anisotropy energy, the single one the conduction electrons, and the solid circles stand 
for the exchange interaction. 



FIG. 6. The second order self-energy correction for the impurity spin. The double line repre- 
sents the impurity spin with the anisotropy energy, the single one the conduction electrons, and 
the solid circles stand for the exchange interaction. 
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FIG. 9. The resistivity for (a) 5 = 2 and (b) S = 5/2 for different values of K. (1) K = 
K = 0.02 K (3) = 0.05 K (4) K = 0.1 K (5) K = 0.2 K {Q) K = 0.5 K (7) = 1 K (8) 
K. The initial parameters were choosen as jq = 0.0435 and = 10^ K, Tk = 0.3 K. 
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FIG. 10. The resistivity for (a) 5 = 2 and (b) S = 5/2 for different values of t/a. (1) 



t/a = oo{K = 0) (2) t/a = 200^ (3) t/a = 100^ (4) t/a = 50^ (5) t/a = 25^ (6) t/a = 10 



K 



(7) t/a = 6i. The initial parameters were choosen as jo = 0.0435 and Dq = 10^ K, Tk = 0.3 K. 
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FIG. 11. The calculated coefHtient B as a function oit/a. The Kondo temperature was choosen 
as Tk = 0.3 K. 




1000 2000 „ 3000 4000 5000 

t(A) 

FIG. 12. Fit on the experimental data (squares) by the calculated formula B{t) = p^^^ Bcaidt / a) 
(the Kondo temperature was choosen as Tk = 0.3 K) with fitting parameters p^^^ = 20nficm, 
a = 247. 7A K (solid line) and by the phenomenological theory B(t) = i?oo(l — 2A/i) with fitting 
parameters Boo = 4.87nf2cm, A = 207.5A (dashed fine). The fit is not too sensitive to small 
changes (< 5%) in a. 
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FIG. 13. The effect of the mean free path on the Kondo resistivity in the presence of anisotropy 
arising only one of the surface in a thin film with thickness L = 600 A for5' = 2atT = 0.6K (1) 
K = {2) K = A^e'"^/^. The Kondo temperature was choosen as Tk = 0.3 K. 
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FIG. 14. Bilayer structure 



FIG. 15. Stripes with same thickness t and changing width w. 
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